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ABSTRACT 

The bulk Lorentz factor of the gamma-ray burst (GRB) cjecta (Fq) is a key parameter to understand 
the GRB physics. Liang et aL have discovered a correlation between Fq and isotropic 7-ray energy: 
Fo oc -E°fso52- By including more GRBs with updated data and more methods to derive Fq, we 



confirm this correlation and obtain Fq 



01 fO.29 
^-•^ -^7,130, 52- 



Evaluating the mean isotropic 7-ray luminosities 



of the GRBs in the same sample, we discover an even tighter correlation Fg 



249L0;30^_52- We 

propose an interpretation to this later correlation. Invoking a neutrino-cooled hyperaccretion disk 
around a stellar mass black hole as the central engine of GRBs, we derive jet luminosity powered by 
neutrino annihilation and baryon loading from a neutrino-driven wind. Applying beaming correction, 
we finally derive Fq oc i^ fg^^, which is consistent with the data. This suggests that the central engine 
of long GRBs is likely a stellar mass black hole surrounded by a hyper-accreting disk. 

Subject headings: gamma-ray: bursts, Lorentz factors - accretion, accretion disks -black hole physics 
- neutrinos 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) are among the r nost pow- 
erful explosions in the universe (jPiranl I20Q4I : iMeszarosI 
120061: lZhan^l2007[ ). It is well known that GRBs are pro- 
duced by relativistic outflows. The bulk Lorentz factor 
during the prompt GRB emission phase (Fq, also called 
"initial" Lorentz factor to be differentiated from the de- 
caying Lorentz factor during the afterglow phase) is a 
very important parameter to understand the physics of 
GRBs. There have been several methods to infer Fq: (1) 
Taking the peak time of the early afterglow light curve 
as the deceleration time of the external forward shock, 
one can estimate Fq, which is twice of the Lorentz factor 
at the deceleration time (Sari & Piran 1999). A non- 
detection of such an afterglow peak time due to a late 
response time or contamination of other e mission com- 
pone nts can lead to a lower limit of Fq (e.g. IZhang et al. 
2006!) . (2) The "compactness problem" constraint (jPiran 
1999), i.e., the requirement that GRBs are optically 
thin to two photo n pair product i on also yiel ds a lower 
limit on Fn (Lithwick fc Sari| [200ll: iGupta fc Zhang.„2008i : 
lAbdo et all I2009allbl lcl). (3) During the prompt emis- 
sion phase, the ex ternal shock is already growing (e.g. 
IMaxham fc ZhanEj r2009). An upper limit of Fq can be 
derived from the data based on the requirement that the 
external shock emission is n ot bright enoug h during the 
prompt emission phase (Zo u fc Piranll20fOI ). 

Constraining Fq of GRBs and studying their statisti- 



^ School of Physics, Huazhong University of Science and Tech- 
nology, Wuhan, 430074, China. Email: zouyc@hust.edu.cn 
(YCZ), leiwh@hust.edu.cn (WHL) 

^ Department of Physics and Astronomy, University of Nevada 
Las Vegas, 4505 Maryland Parkway, Box 454002, Las Vegas, NV 
89154-4002, USA. Email: zhang@physics.unlv.edu (BZ) 

^ Department of Physics, Guangxi University, Nanning 
530004, China 

* GXU-NAOC Center for Astrophysics and Space Sciences, 
Nanning, Guangxi 530004, China 



cal properties is essential to constrain the ph ysical origin 
of GRB prompt emission ([Liang et al.ll2010f ). In particu- 
lar, different theoretical models demand different correla- 
tions between Fq and E^^iso or i^.iso (Zhang & Mcszaroi 
l2002f ) in or der to account for the observed Ep — i?^,iso 
correlation (jAmati et al.l 120021 ) . By constraining Fq of 
about 20 GRBs that show the deceleration featu re in 
the early afterglow lightcurves, ILiang et"aLl (|2010f ) dis- 
covered a tight correlation between Fq and Ej^iso, i-e. 
Fo ~ 182(£:^.iso/10^^er g)°-^^ Using a different method 
to derive Fq, .Ghirlanda et al.l ()2011[ ) confirmed a corre- 
lation between Fg and E^^iso, but with a different power 
inde^l^- 

In this paper, we work on an expanded sample and ap- 
ply more methods to constrain Fq for about 50 GRBs. 
W e test and confirm t he To — E^^iso correlation discovered 
bv ILiang et al.l (|2010D . and also investigate a Fq — L^.iso 
correlation, where L-y^so is the mean luminosity of the 
burst. Since the Fq — E^^iso c orrelation was no t in- 
terpreted in the previous work IjLiang et al.l 120101) , we 
also attempt to propose an interpretation in this pa- 
per. We derive jet luminosity and baryon loading from 
a black hole - neutrino-cooling-dominated-flow (NDAF) 
disk central engine model, and find that this central en- 
gine model can naturally account for the Fq — i-y,iso cor- 
relation discovered in this paper, an d hence, can also in - 
terpret the Fq — i^-y, iso correlation of ILiang et al.l (|2010f) . 

We arrange this paper as follows. In Section 2, the 
methods of Fq derivations based on three methods are 

^ The Ghirlanda et al. (2011) method applies the Blandford- 
Mckee (BM) self-similar deceleration solution (Bla ndford McKe^ 
119761) and extrapolates it backwards to derive Fq- However, around 
the deceleration stage, the dynamics has not entered the BM self- 
semilar solution yet. Also the intersection of the two asymptotic 
power law phases (as adopted by Ghirlanda et al. 2011) may not 
correspond to the observed peak time of afterglow light curve. We 
regard the Ghirlanda et al. (2011) method not more precise than 
the conventional method, and still adopt the conventional method 
to derive Fq in this paper. 
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summarized. We then apply the methods to the available 
GRBs on which these methods can be used, and present 
the To — E^^iso and Tq — L^^iso correlations in Section 3. 
In Section 4, a physical interpretation to the Tq — L^^iso 
correlation is presented. Our results are summarized in 
Section 5 with some discussion. 

2. METHODS OF CONSTRAINING To 

We apply three methods to co nstrain F p, namely, 

(A) the afterglow onset method (iSari fc Pir an 1999), 

(B) pair opacity constraint method ( Lithwick &: Saril 



and 



l200ip, and (C) ear ly external forward emission method 
(IZou fc Piranl[20Tl . 

Method A is the most common method, which uses the 
peak of the early afterglow light curve to determine the 
deceleration time of the external forward shock. In the 
so-called "thin shell" regime, the initial Lorentz factor 
Tq is twice of the Lorentz factor at the deceleration time. 
For a constant density medium, one has 



1.4 



3^^,iso(l + 

327rnmpc577i3^^i^ 



-1 1/8 



(1) 



where n is the medium number density, rrip is the proton 
rest mass, 77 is the ratio between the isotropic gamma- 
ray energy and the isotropic blast wave kinetic energy, 
and tpeak is peak time of the afterglow, which is also 
taken as the deceleration time. The derived Fq is rather 
insensitive to n and 77, but mildly depends on tpeak (-3/8 
power). If the peak time is not detected, ipoak is regarded 
to be prior to the earliest afterglow observing time (i.e., 
iobs > ^poak)- This givcs a lower limit on Fq. Notice 
that in Eq.(IT]), we have taken the coefficient as 1.4 rather 
than the commonly used 2. This more precise factor 
comes from two factors: First, the deceleration radius 
is defined by the condition M = Afo/Fdoc rather than 
M = Mq/Tq (where M is the shocked ISM mass, and Mq 
is the original mass of the ejecta), since at this radius, 
the shocked ISM and the ejecta have the same inertia. 
Second, instead of adopting rdoc — 2F^^^cidoc, we apply a 
differential form dr ~ 2T^cdt, and numerically integrate 
it from t — to t = t^cc to get rdoc = 4.4F^^j.cidcc- Here 
rdec is the deceleration radius and tdcc is the deceleration 
time, which also corresponds to the peak time tpcak- 

Method B requires that observed high energy 7- 
rays (e.g. those in the GeV range) are optically thin 
to electron-positron pair production with softer tar- 
get photons in the emission region. This yields a 
lower limit on the Lorentz factor of the emitting re- 
gion (jLithwick fc Saril l200l[ l^. The lower limit can be 
obtained by requiring that the observed highest energy 
photons with energy i?max have an optical depth smaller 
than unity: 



Fo > f 2(/3 + 2) 



(1 + z)- 



(2) 



® This method makes the assumption that the GRB emission 
radius -Rgrb is related to Fq via -Rgrb — Fgc5T. Some GRB 
prompt emis sion models do not satisfy such a condition (e.g. 
INaravan fcK umar 2009; Zhang & Yan 201^). The lower limit of 
Fo cannot be uniquely de rive d, since the cutoff energy is a function 
of both Fo and Rgrb (Gupta fc Z hang 20M)- 



f = 2.1 X 10' 



(D/7Gpc)2(0.511)-'3+Vi 



(<5r/0.1s)(/3-l) 



(3) 



where f3 is the photon spectral index in the MeV band, 
with a typical value between 2 and 3, D is the luminos- 
ity distance, ST is the minimum variability time scale of 
the prompt emission, and /i is the observed number of 
photo ns per second per cm^ per MeV at the energy of 1 
MeV (jLithwick fc Saril |2001[) . We notice that there are 
a few bursts whose Fq's constrained using this method 
are inconsistent with those derived from other two meth- 
ods. Instead, we apply a modified version of Method B, 
which assumes the high energy emission and the prompt 
MeV emission are from two diffe rent emitting regions 
(|Zou et al.ll2"0Tlt IZhao et al.ll20Tll) . 

Method C considers the quiescent periods between the 
prompt emission pulses, in which the signal of external 
shock has to go down the instrume nt thresholds. Thi s 
would place an upper limit on Fq ( Zou fc PiranI [20T0I) . 
The constraint of Fq for a uniform density medium is 



Fo < 340(l + z)3/,,lim,-28«i?i 



,-l''20 



-V(i. 



where /i/jim ~ 10~^^erg cm^^s^^Hz^^ is the limiting flux 
density of the observing instrument, Y is the Compton 
parameter for synchrotron self- Compton scattering, ec 
is the equipartition factor for internal energy density of 
electrons, eb is the equipartition factor for the magnetic 
energy density, and is the first quiescent time in the 
observer's frame, and v is the observing frequency. In 
this paper, we take the conventional notation Q = Qk x 
10*'' if not specified. 

3. SAMPLE SELEGTION AND CORRELATIONS 

Using the methods above, we can constrain Fq for the 
bursts with enough observational data. The parame- 
ters of 51 GRBs in our sample are presented in Table 1, 
which include spectroscopically confirmed redshift (z), 
burst duration (Tgo), derived initial Lorentz factor Tq, 
isotropic 7-ray energy (Ey^iso), and isotropic mean 7-ray 
luminosity {L.yAso = (1 + z)Ey^iso/TQo). Within the sam- 
ple, 38 GRBs have Fq calculated using Method A (Refs 
a, b and d in Table 1). As methods B and C can only 
get a range for the derived Lorentz factor, the fit for the 
relations of Fq — E^ iso and Fq — iso are from these 38 
GRBs only. 

With the data listed in Table 1, a correlation analysis 
between logFo and log L^, iso data set yields a Pearsons 
correlation coefficient with C = 0.79, which is tighter 
than the logFo — log i?^, iso correlation with ( — 0.67. 
We plot Fq versus i?-y,iso and L^^iso in Fig. 1 and Fig. 2, 
respectively. Visibly one can see a strong correlation in 
both plots. The best fitting results are: 

logFo = (1.96 ±0.002) + (0.29 ±0.002) log £:.^,iso,52 (5) 

with ( = 0.67, and 

logFo = (2.40 ± 0.002) + (0.30 ± 0.002) logL.^,iso,52 (6) 

with C = 0.79. 
These correlations can be translated to 
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and 



(8) 



I t can be seen that t he Fq— _B-yjso correlation discovered 
by ILiang et al.l (j201 (y) is confirmed. The smaller coeffi- 
cient (91 instead of 182) is mainly caused by the smaller 
(but more precise) factor 1.4 (rather than 2) in Eq.([T]). 
We also found a tighter Fq — i-y.iso correlation, suggest- 
ing that it may be more intrinsic than the Fq — E^/iso 
correlation. As Fq, i?7,iso and L^y^iso are all z-dependent 
quantities, there might be a selection effect involve d so 
that the correlation may not be intrinsic (fButler et al.l 
|2009( ). In order to test this possibility, we study the 
Fq — L-y^iso relation with the following procedure: 1. 
We randomly produce a set of r edshifts according to 
the G RB z-distribution given by iWanderman fc PiranI 
(|2OT0l) : 2. assign these random artificial redshifts to the 
bursts to replace the observed ones; 3. calculate the Fq 
and L^^iso according to the artificial redshifts; 4. cal- 
culate the correlation coefficient ( of logFo — logL^y^so 
correlation for each realization; 5. redo step 1 through 
4 10000 times, and get a distribution of correlation co- 
efficient; 6. compare the most probable coefficient with 
the coefficient generated from the real data. The most 
probable coefficient from our simulations is 0.63, which is 
clearly smaller than the one derived from the real data, 
C = 0.79. This means that the Fq — L^^iso relation is 
likely intrinsic, not caused by a selection effect from z- 
dependence parameters. 

We notice two outliers to both correlations: 
GRB060614 and GRB080129, whose Fq's are de- 
rived using Method A from a late optical bump, which 
lead to Fq < 100 for both cases. It is possible that these 
bumps a re cau sed by other mechanisms (e.g. energy 
injection, iXu et al. 2009). If this is the case, the derived 
Fo for the two bursts can be regarded as lower limits. 

4. THEORETICAL INTERPRETATION 

The most popular model of GRB central engine in- 
vokes a stellar mass black hole su rrounded by a hyper- 
accre t ing disk (e.g. IPopham et a l. 1999; Narayan ct al. 
20011: iDi M atteo et al.l l20 02t iKohri fc Miiieshige .2002^ 
Gu et al.ll20 06: Chen & Bc loborodovll2007l: iJaniuk et al.l 
2007tlLei et al.il2009ll . In the inner region of such a hy- 
peraccretion disk a large amount of energetic neutrinos 
are emitted, carrying away the viscous dissipation en- 
ergy of the accreted gas. If the accretion rate is not too 
low, neutrino annihilation {vi' e^e~) can launch a 
relativistic jet powerful enough to account for the GRB. 

For a system with black hole mass M and spin a* , the 
neutrino annihilation power Ei^p from the hyperaccretion 
disk depends on the accretion rate M (for Mign < M < 
Mtrap) as (jZalamea fc Beloborodovll2011f ). 



1.1 X W^x- 



m 



erg s 



(9) 



where M3 = M/3Mq, m = M/MqS'^, x^s = 
rms{a*)/rg, and Vg — 2GM/c^. Here rms is the radius 
of the marginally stable orbit, wh ich is a function of 
the black hole spin a* (Page fc Thornelll974[ ). We have 
Xras — 0.97 for a* — 0.95. The t wo critical accretion rates 
Mign and Mtrap are defined in iZalamea fc Beloborodovl 



high enough to ignite neutrino emitting reactions. If 
M > A'/trap , the emitted neutrinos become trapped in 
the disc and advected into the black hole. For the disk 
with viscosity a = 0.1, we find Mjgn = 0.071AfQS~^ and 

Mt.^p = 9.3M0S-1 for a* = 0, and Mig„ = O.O2IM0S-1 

and Mtrap = I.SAIqS^^ for a* ~ 0.95. 

Most neutrino annihilation energy is converted into ki- 
netic energy of baryons after acceleration, and the jet 
reaches a Lorentz factor 



M„c2 



(10) 



where Mi, is the neutrino-driven mass loss rate from the 
disk. The mass loss rate M^ is related to the total neu- 
trino power E^ through (jMetzger et al.l [20081 ) 

M. ^ 10-'El%{el,f/\l^'M^'ih/r)-'MQS-' (11) 

where re = r/lO'^cm, £1,^52 = i^iz/lO^^erg s~^, e^, = eio x 
lOMeV is the mean energy of neutrinos, and h is the half- 
thickness of disk. For a^, — 0.95, the total neutrino power 
from the disk is Ejy ~ 0.15Mc2 (jChen fc Beloborodovl 
[20f)l . 

For a neutrino dominated accretion flow (NDAF) , both 
(which is a function of disk temperature) and h are 
independent of the accretion rate rh. This result can be 
checked with the analytical solution o f hyper-accreting 
disk obtained bv IPopham et al.l (|1999D (i.e., their equa- 
tions (5.3)and (5.4)). So, based on Eq. ([TT1) . the mass loss 

rate M^, is just related to m as cx nr'^. Com- 

bining this dependence with Eq.Q, one drives M^, (x 
m^/^ oc E^J'^'^ . And then inserting it to Eq. lTUl) . we 

therefore obtains Fq oc Ey^, / Mi, oc Eji, . 

The relativistic jet with Lorentz factor Fq will dissipate 
its kinetic energy via internal shocks with efficiency 77 and 
produce gamma-ray emission, i.e., ~ "q^E^fj. Assum- 



7/27 



In 



ing a constant rj^ for all GRBs, one can get F oc L 
order to connect and L-yjso, one needs to further take 
into account the beaming correction, i.e. = fbL^,iso- 



One then gets 



fb ^7 fb ViEuv oc /f, E^ 



where <C 1 is the beaming factor. 

The general dependence of fb on the properties of cen- 
tral engine is unknown. However, one can gain insight 
directly from observations. Following Amati et al. (2002, 
2006, 2008), the relationship between the isotropic equiv- 
alent energy radiated during the prompt phase (i?-y,iso) 
and the rest-frame peak energy in the GRB spectrum 
{E'p ) is E'p (X i?°;fso- By combining it with the Girlanda 
relation E^ oc (E'p)^/^ (Ghirlanda et al. 2004), where the 
beaming-corrected energy Ej = /fc-E^jso, we obtain the 
relation between fb and Ej^iso as fb oc E^'^^l^^. Since 
L^.iso oc Ej^iso, we get fb oc L'^^sl^^- One can see that fb 
is very insensitive to L^^iso and i^^^iso- 

Now we can obtain the relation between Lorentz factor 
Fq and the isotropic luminosity L-yjso based on the above 
scalings, i.e. 



Fq oc eH^'^ oc (fbL^ 



-.7/27 



(12) 



(|2011f) . If M < Mign, the disc temperature is not In view of the large scatter of the applied empirical 
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Amati- and Ghirlanda-correlations, we regard that this 
theoretically motivated correlation agrees with the sta- 
tistical correlation ([8]). 

5. CONCLUSIONS AND DISCUSSION 

By including more recent GRBs and by engaging 
more methods to constrain Fq, we have critically re- 
analyzed the statis tical correlation between Fq and Eiso 
(jLiang et al.ll201Cll ). We confirmed the correlation and 
found Fq ~ 91ii^° ?Q Furthermore, we found an even 
tighter correlation between Fq and the mean isotropic 
7-ray luminosity, which reads Fq — 249L'^ fJ^ i^2- 

We also proposed an interpretation to the Fq ^ L^so''^ 
correlation within the framework of a black hole - NDAF 
disk GRB central engine model. By invoking a neutrino- 
annihilation powered jet and by calculating baryon load- 
ing from a neutrino-driven wind, we get a Fq cx: L^^^ 
correlation. Further considering the beaming factor /f,, 
which is insensitive to L^^so as evidenced from the em- 
pirical Amati and Ghirlanda correlations, we finally de- 
rived Fq ^ i^'.iso- In view of the large scatter of various 
correlation, we regard that this model prediction is well 
consistent with the observed Fq — i-yjso correlation. 

The existence of the Fq — L^^iso and Fq — iJ^.iso cor- 
relations and the success of interpreting them within 
the black hole - NDAF central engine model hint that 
the GRB central engine is likely a hyper-accreting black 
hole. The interpretation invokes a neutrino-annihilation- 
powered jet, which is justified for a reasonably high accre- 
tion rate and a not very rapid black hole spin (W.-H. Lei 
& B. Zhang 2011, in preparation). On the other hand, 
recently arguments have been raised to support a mag- 
netically dominated jet from GRB s (e.g. iZhang fc Pe'ed 
[20091: IFmI[2010I: IZhang fc Ya^[20Tll) . Studies of the black 



hole central engine models also s uggest that ma gnetic 
fields play an important role (e.g. iLei et al.|[2009l ). The 
baryon loading process in a magnetically dominated jet is 
more complicated, and has not been studied carefully in 
the literature. Whether the Fq — L^^iso correlation can be 
still interpreted in a magnetized black hole cen tral engine 
model (e.g. Blandford & Znaick 1977; Mcsza ros fc Reed 
119971: iWaiig et al.,,2002; .Yuan fc Zhana.2011.) is subject 
to further investigations. 

Since only one short GRB (090510) is included in our 
sample, our correlations and interpretation apply to long 
GRBs only. 

Recently, iWu et all ()201lD discovered an intriguing 
universal correlation between synchrotron luminosity 
and Doppler factor for GRBs and blazars. Our interpre- 
tation cannot be extended to blazars, since the accretion 
rate inferred from blazars are not in the NDAF regime. If 
indeed the two phenomenon share the same physics, then 
the correlation may stem from a more profound physical 
origin, which is beyond the scope of this paper. 
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Fig. 1. — The plot of the derived initial Lorentz factor Fq vs. the isotropic equivalent 7-ray energy -E-^,iso. The solid line, Fq ~ ^^E^'f^a 52' 
is the best fit to the derived values, which are in solid circles. The Pearson's correlation coefficient is 1^ = 0.67. The triangles are the bursts 
having only lower limits and the range segment are the bursts having upper and lower limits listed in Table 1, which are not included in 
the fitting process. The star is the only short burst GRB 090510. 
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Fig. 2. — The initial Lorentz factor Fq vs. the isotropic equivalent 'y-ray luminosity L^ ^so- The best fitting line is Tq ~ 2491°-?°^ 52 with 
correlation coefficient = 0.79. The notations are the same as in Fig. 1. 
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TABLE 1 

The quantities of the GRBs in our sample. Main references are listed in the last 
COLUMN. Others are marked directly after the quantities. 
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i'o 
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1.6lli>^l 
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2.3 ± 0.52(441 
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39.18(50] 


0.71 ±0.071 


d 


100621A 


0.542(43] 


52 


4.37 ±0.5 


63.6 ±1.7(73] 


0.11 ±0.015 


d 


100728B 


2.106(43] 


373 


3 ±0.3 


12.1 ±2.4(74] 


0.77 ±0.23 


d 


100906A 


1.727(43] 


369 


33.4 ±3 


114.4 ±1.6(75] 


0.8 ±0.083 


d 


110205A 


2.22(43] 


177 


56 ±6 


257 ± 25(76] 


0.7 ±0.14 


d 


110213A 


1,46143] 


223 


6.4 ±0.6 


48 ± 16(77] 


0.33 ± 0.14 


d 



Ref eren ces. — [1] lAkerlof eFaLl l[l999l ): [2] IGalama et"an CS^; [3 ] IVreeswiik et al.l 

HWOB): [4 1 IWicrscma ct aE ('2005): [5] Dc iPasqualc ct al." (2006); [61 FolcY Ct al. (2005' ); |7| 
IRol et al.l (|2005); [81 Dc Pasqualc ct al. (2007); [91 Lcdoux ct al. (2005); [101 Jakobsson e taU 

(120051): [111 Cucchiara ct al. (2006); [121 Molinari ct al. (2007); [13] Falcone ct al. (2006); [141 
I Ferrero et a l. (2009); [15] Page ct al. (2007); [16] MundclK 2007); [17] Fugazza ct al. (2006); [18] 
IRoI et al.' (^2006): [19] Jakobsson ct al. (2007a); [20] Bloom (2006); [21] Jaunscn ct al. (2007); 

g2l Chen ct al. (2007); [23] Jakobsson ct al. (2007a); [24] Ccnko ct al. (2007b); [25] Pcrlcv ct al] 

il2008al); [26 ] Prochaska ct al. (2007); [27] Ccnko ct al. (2007a); [28] Lcdoux ct al. (2007); [29] 
IGreiner etal . (2009); [30] Stratta ct al. (2009); [31] Ran et al. (2009); [32] Wierscma etall 

120081); |33T Cucchiara (2008); [34] Fynbo ct al. (2009); [35] Guidorzi et al. (2009b); [361 
[Krlcv et al. (2008); [37] Prochaska ct al. (2008); [38] Abdo ct al. (2009a); [39] Landsman etahl 

(poos'!: [401 .Mela ndri ct al. (2010); [41] Ccnko et al. (2011); [42] Wierscma et al. (2005) ; [431 
IGhirlanda et al.l ( 2011); [44] Crew ct al. (2003); [45] Donaehy ct al. (2006); [46] Blustin etaU 

(|2006l): I47| |Pandey ct al. (2006); [48] Rykoff et al. (2006); [49] Ccnko ct al. (2006); [501 
ISakamoto et al. (2011); [51] Curran ct al. (2007); [52] Ziaccpour et al. (2008); [53] Gchrcls ct at] 

(120061): 154 ] Klotz ct al. (2008); [55] Coyino ct al. (2010); [56] Schady ct al. (2007); [57f 
ITroia et all (2007): [58] Chester ct al. (2008); [59] Ferrcro ct al. (2008); [60] Melandri et all 

(120091): 161 ] Kong ct al. (2010); [62] Wang et al. (2008); [63] Racusin ct al. (2008); [64] 
lYuan et all (=2008): [65] Martin-Carrillo ct al. (2008); [66] Kriihler ct al. (2009); [67] Page et^ 

(120091): 1681 dc Ugartc Postigo ct al. (2010); [69] Dc Pasqualc et al. (2010); [70r Abdo ct al] 
JgOOgd); [711 Swenson ct al. (2010); [72] Grubcr ct al. (2011); [73] Ukwatta ct al. (2010); [74] 
[Barthclmy ct al. (2010 a); [ 75] Barthclmy c t al. (2 010b ); [76] Markwardt ct al (2011); [77] 

IRni-tViolm^r et n1 I (OTVUtl- fa\ \Mc\anr\r\ ot nl I IWlTiTl- (h\ IT.inno- et al I /f9m71^• (n\ I7nii .lir PirsjTil 



